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The problem of counting tilings of a plane region using specified tiles can often be recast as the problem of 
counting (perfect) matchings of some subgraph of an Aztec diamond graph A,,, or more generally calculating 
the sum of the weights of all the matchings, where the weight of a matching is equal to the product of the 
(pre-assigned) weights of the constituent edges (assumed to be non-negative). This article presents efficient 
algorithms that work in this context to solve three problems: finding the sum of the weights of the matchings 
of a weighted Aztec diamond graph A„; computing the probability that a randomly-chosen matching of A„ will 
include a particular edge (where the probability of a matching is proportional to its weight); and generating 
a matching of A„ at random. The first of these algorithms is equivalent to a special case of Mihai Ciucu's 
cellular complementation algorithm |^ and can be used to solve many of the same problems. The second of 
the three algorithms is a generalization of not-yet-published work of Alexandru lonescu, and can be employed 
to prove an identity governing a three- variable generating function whose coefficients are all the edge-inclusion 
probabilities; this formula has been used as the basis for asymptotic formulas for these probabilities, but a 
proof of the generating function identity has not hitherto been published. The third of the three algorithms is a 
generalization of the domino-shuffling algorithm presented in [0]; it enables one to generate random "diabolo- 
tilings of fortresses" and thereby to make intriguing inferences about their asymptotic behavior. 



1. INTRODUCTION 



matching of an Aztec diamond graph of order 3: 



1.1. Background 



Let L be the rotated square lattice with vertex set : 
ij e Z, ! + i is odd}, and with an edge joining two vertices 
in the graph if and only if the Euclidean distance between the 
corresponding points in the plane is ^/2. Define the Aztec 
diamond graph of order n (denoted by A„) as the induced 
subgraph of L with vertices satisfying — n < i,j < n, as 
shown below for n = 3. 




A partial matching of a graph is a set of vertex-disjoint edges 
belonging to the graph, and a perfect matching is a partial 
matching with the property that every vertex of the graph be- 
longs to exactly one edge in the matching. Here is a perfect 



(Hereafter, the term "matching", used without a modifier, will 
denote a perfect matching.) 

Matchings of Aztec diamond graphs were studied in in 
the dual guise of domino-tilings (see below for the domino- 
tiling dual to the previously depicted matching). 
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Note that this is rotated by 45 degrees from the original pic- 
ture presented in That article showed, by means of four 
different proofs, that the number of matchings of an Aztec 
diamond graph of order n is 2"'"+''/^ (as had been conjec- 
tured a decade earlier in the physics literature [[l0[]). One 
of the proofs of this result (developed by authors Kuper- 
berg and Propp) used a geometric-combinatorial procedure 
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dubbed "shuffling" (short for "domino-shuffling"). Shuffling 
was originally introduced purely for purposes of counting the 
domino-tilings of the Aztec diamond of order n, but it later 
turned out to be useful as the basis for an algorithm for sam- 
pling from the uniform distribution on this set of tilings. An 
undergraduate, Sameera Iyengar, was the first person to imple- 
ment domino-shuffling on a computer, and the output of her 
program suggested that domino-tilings of Aztec diamonds ex- 
hibit a spatially-expressed phase transition, where the bound- 
ary between "frozen" and "non-frozen" regions in a random 
domino-tiling of a large Aztec diamond is roughly circular 
in shape. A rigorous analysis of the behavior of the shuf- 
fling algorithm led to the first proof (by Jockusch, Propp, and 
Shor 1 1 1 1) of the asymptotic circularity of the boundary of the 
frozen region (the "arctic circle theorem"). 

Meanwhile, Alexandru lonescu, also an undergraduate at 
the time, discovered a recurrence relation related to shuffling 
that permits one to calculate fairly efficiently, for any edge e 
in an (unweighted) Aztec diamond graph, the probability that 
a random matching of the graph contains the edge e. This 
theorem allowed Gessel, lonescu, and Propp (in unpublished 
work) to prove a conjecture of Jockusch concerning the value 
of this probability when e is one of the four central edges in the 
graph. More significantly, the theorem allowed Cohn, Elkies, 
and Propp ||^] to give a detailed analysis of the asymptotic 
behavior of this probability as a function of the position of the 
edge e as the size of the Aztec diamond graph goes to infinity. 

During this same period, it had become clear that 
(weighted) enumeration of matchings of weighted Aztec di- 
amond graphs had relevance to tiling-models other than 
domino-tiling. Notably, Bo- Yin Yang had shown [ p2| ] that 
the number of "diabolo-tilings" of a "fortress of order n" was 
given by a certain formula conjectured by Propp (for more on 
this, see section 8). Yang's proof made use of Kuperberg's re- 
formulation of the problem as a question concerning weighted 
enumeration of matchings of a weighted version of the Aztec 
diamond graph in which some edges had weight 1 while other 
edges had weight j (see the end of section 5 of this article 
and the longer discussion in section 8). It was clearly desir- 
able to try to extend the shuffling algorithm to the context of 
weighted matchings, or as physicists would call it, the dimer 
model in the presence of non-uniform bond-weights, but at the 
time it was unclear how to devise the right extension. Ciucu's 
work iQl went part of the way towards this, and in particular it 
gave a very clean combinatorial proof of Yang's result. How- 
ever, it was not clear whether one could efficiently generate 
diabolo-tilings using some extension of domino-shuffling. 

The purpose of the present article is to give the details of 
just such a generalization of domino-shuffling into the general 
context of Aztec diamond graphs with weighted edges. At the 
same time, my goal is to remove some of the mystery that 
hung about the shuffling algorithm in the original exposition 
in [|]|, and to give a exposition of the work of lonescu, which 
up until now has not been described in the literature at any 
useful level of detail. It should be stressed 

• that the counting algorithm (presented in section 2 and 
analyzed in section 5) is essentially a restatement of Mi- 
hai Ciucu's cellular graph complementation algorithm 



• that the probability computation algorithm (presented 
in section 3 and analyzed in section 6) is a reasonably 
straightforward generalization of lonescu's unpublished 
algorithm for the unweighted case, 

• and that the random generation algorithm (presented in 
section 4 and analyzed in section 7), although alge- 
braically more complicated, represents no conceptual 
advance over the unweighted version of shuffling pre- 
sented by Kuperberg and Propp. 

I think that the major contribution offered by this paper is not 
any one of these algorithms in isolation but the way they fit 
together in a uniform framework. 

The essential tool that makes generalized shuffling work is 
a principle first communicated to me by Kuperberg, which 
is very similar to tricks common in the literature on exactly 
solvable lattice models. This principle is a lemma that as- 
serts a relationship between the dimer model on one graph 
and the dimer model on another, where the two graphs differ 
only in that a small patch of one graph (sometimes called a 
"city") has been slightly modified in a particular way (a pro- 
cess sometimes called "urban renewal"). Urban renewal is a 
powerful trick, especially when used in combination with an 
even more trivial trick called vertex splitting/merging. (Such 
graph-rewriting rules are not new to the subject of enumer- 
ation of matchings; see for instance the article of Colbourn, 
Provan, and Vertigan for a very general approach that 
uses the "wye-delta" lemma in place of the "urban renewal" 
lemma.) Urban renewal, applied in locations in an Aztec 
diamond of order n, essentially converts it into an Aztec di- 
amond of order n — 1 . In this way, a relation is established 
between matchings of A„ and matchings of j . In partic- 
ular, one can reduce the weighted enumeration of matchings 
of a weighted Aztec diamond graph of order n to a (differ- 
ently) weighted enumeration of matchings of the Aztec dia- 
mond graph of order n — 1 . 

The results described in this article were first commu- 
nicated by means of several long email messages I sent 
out to a number of colleagues in 1996. I am greatly 
indebted to Greg Kuperberg, who in the Fall of 2000 
turned these messages into the first draft of this article, 
and who in particular created most of the illustrations 
(thereby introducing me to the joys of pst ricks). Also, 
Harald Helfgott's program ren (currently available from 
http : / / www . math . wise . edu/ ^propp/ ; see the files 



ren . c, ren . h, and ren . html) was the first implementa- 
tion of the general algorithm, and it facilitated several useful 
discoveries, including phenomena associated with fortresses 
||l4|]. Helfgott and lonescu and Iyengar were three among 
many undergraduate research assistants whose participation 
has contributed to my work on tilings; the other members 
of the Tilings Research Group (most of whom were under- 
graduates at the time the research was done) were Pramod 
Achar, Karen Acquista, Josie Ammer, Federico Ardila, Rob 
Blau, Matt Blum, Ruth Britto-Pacumio, Constantin Chiscanu, 
Henry Cohn, Chris Douglas, Edward Early, Marisa Gioioso, 
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David Gupta, Sharon Hollander, Julia Khodor, Neelakan- 
tan Krishnaswami, Eric Kuo, Ching Law, Andrew Menard, 
Alyce Moy, Ben Raphael, Vis Taraz, Jordan Weitz, Ben 
Wieland, Lauren Williams, David Wilson, Jessica Wong, Ja- 
son Woolever, and Laurence Yogman. I thank Chen Zeng, 
who has made use of the generalized shuffling algorithm [ p3[ | 
in his own work on and thereby encouraged me to write up this 
algorithm for publication. Lastly, I thank the anonymous ref- 
eree, whose careful reading and helpful comments have made 
this a better article in every way. 



1.2. Two examples 

Here are two illustrations of the flexibility of the procedures 
described in this article. The figure below shows a subgraph 
of the Aztec diamond graph Aj. This can be thought of as 
an edge-weighting of the Aztec diamond graph of order 5 in 
which the marked edges are assigned weight 1 and the remain- 
ing (absent) edges are assigned weight 0. 



c = 2. 




Let V be the set of triangles, and let E be the set of pairs of 
triangles sharing an edge. (Pictorially, we may represent the 
elements of V by dots centered in the respective triangles, and 
the elements of E by segments joining elements of V at unit 
distance.) The graph {V,E) is the hexagon honeycomb graph 
b c- Here, for instance, is the hexagon honeycomb graph 

^2,2,2- 




It is easy to see that the twelve isolated edges all are "forced", 
in the sense that every matching of the graph must contain 
them. When these forced edges are pruned from the graph, 
what remains is the 6-by-6 square grid. More generally, for 
any n, one can assign a 0, 1 -weighting to the edges of the 
Aztec diamond graph of order 2n — 1 so that the matchings 
of positive weight (i.e., weight 1) correspond to the match- 
ings of the 2«-by-2« square grid. One can use this correspon- 
dence in order to count matchings of the square grid, deter- 
mine inclusion-probabilities for individual edges, and gener- 
ate random matchings (for a technical caveat, see section 9, 
item 1). 

The second example concerns matchings of graphs called 
"hexagon honeycomb graphs". These graphs have been stud- 
ied by chemists as examples of generalized benzene-like hy- 
drocarbons 1^, and they also have connections to the study 
of plane partitions [Q]. To construct such a graph, start with 
a hexagon with all internal angles measuring 120 degrees and 
with sides of respective lengths a,b,c,a,b,c, where a, b, and 
c are non-negative integers, and divide it (in the unique way) 
into unit equilateral triangles, as shown below with a = b — 



Perfect matchings of such a graph correspond to tilings of the 
a,b,c,a,b,c hexagon by unit rhombuses with angles of 60 and 
120 degrees, each composed of two of the unit equilateral tri- 
angles. 

The following figure shows, illustratively, how 2 2 
be embedded in the Aztec diamond graph A^: 




More generally, for any a,b,c, one can assign a 0, 1-weighting 
to the edges of a suitably large Aztec diamond graph so that 
the matchings of positive weight (i.e., weight 1) correspond to 
the matchings of H^^ ^ One can use this correspondence in 
order to count matchings of the honeycomb graph, determine 
inclusion-probabilities for individual edges, and generate ran- 
dom matchings (once again subject to the technical caveat dis- 
cussed at the beginning of section 9). 

1.3. Evaluation and comparison of the algorithms 

I have not done a rigorous study of the running-times of 
these algorithms under realistic assumptions, but if one makes 
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the (unrealistic) assumption that arithmetic operations take 
constant time (regardless of the complexity of the numbers 
involved), then it is easy to see that the most straightforward 
implementation of generalized shuffling, when applied in an 
Aztec diamond of order n, takes time roug hly (i gnormg 
factors of logn). We have found that in practice, shuffling is 
extremely efficient. The physicist Zeng has applied the al- 
gorithm to the study of the dimer model in the presence of 
random edge-weights [ ^3| ] . 

Alternative methods of enumerating matchings of graphs 
are legion; the simplest method that applies to all planar 
graphs is the Pfaffian method of Kasteleyn |12|, and there 
is a variant due to Percus JT^ (the permanent-determinant 
method) that applies when the graph is bipartite. For the spe- 
cial case of enumerating matchings of Aztec diamond graphs, 
there are close to a dozen different (or at least different- 
looking) proofs in the literature. One that is particularly wor- 
thy of mention is the recent condensation algorithm of Eric 
Kuo [|l|]. Kuo's algorithm appears to be closely related to the 
algorithm presented here, although I have not worked out the 
details of the correspondence. It is also worth remarking upon 
the resemblance between Kuo's basic bilinear relation and the 
bilinear relation that appears in section 6 of this article; both 
relations are proved in similar ways. 

Once one knows how to count matchings of graphs (or, in 
the weighted case, to sum the weights of all the matchings), 
there is a trivial way to use this to compute edge-inclusion 
probabilities: apply the counting-algorithm to both the orig- 
inal graph and the graph from which a selected edge has 
been deleted (along with its two endpoints). The ratio of 
these (weighted) counts gives the desired probability. If one 
wants to compute only a single edge-inclusion probability, 
this method is quite efficient, and if one uses Kasteleyn's ap- 
proach, one need only calculate two Pfaffians; however, if one 
wants to compute all edge-inclusion probabilities, it is waste- 
ful to compute all the Pfaffians independently of one another, 
since the different subgraphs are all very similar. Wilson [ pi] ] 
shows how one can organize the calculation more efficiently, 
exploiting the similarities between the subgraphs. I believe 
that Wilson's algorithm and the approach given here are likely 
to be roughly equivalent in terms of computational difficulty 
and numerical stability. The new approach has the virtue of 
being much simpler to code. 

Once one knows how to compute edge-inclusion probabili- 
ties, generating a random matching is not hard: one can cycle 
through the set of edges, making decisions about whether to 
include or not include an edge in accordance with the asso- 
ciated inclusion-probability (whose value in general depends 
on the outcome of earlier decisions). The naive way of im- 
plementing this, like the naive way of computing all the in- 
clusion probabilities in parallel, requires many determinant- 
evaluations or Pfaffian-evaluations of highly similar matri- 
ces; as in the case of the previous paragraph, Wilson 
has shown how to use this to increase efficiency, resulting 
in an 0{N^'^) algorithm, where is the number of vertices. 
Note that in the case of Aztec diamonds, = 0{n^) and 
A^^-^ = 0{n^), so this is the same approximate running-time 
as the algorithm given in this article. No algorithm faster than 



Wilson's is known. 



1.4. Overview of article 

The layout of the rest of the article is as follows. Sections 
2 through 7 correspond to the original six installments of the 
1996 e-mail version of the article. Sections 2 through 4 give 
algorithms for weighted enumeration, computation of proba- 
bilities, and random generation; sections 5 through 7 give the 
proofs that these algorithms are valid. The expository strategy 
is to use examples wherever possible, rather than resort to de- 
scriptions of the general situation, especially where this might 
entail cumbersome notation. Section 8 applies the methods of 
the article to the study of diabolo-tilings of fortresses. Sec- 
tion 9 concludes the article with a summary and some open 
problems. 



2. COMPUTING WEIGHT-SUMS: THE ALGORITHM 

Here is the rule (to be immediately followed by an example) 
for finding the sum of the weights of all the matchings of an 
Aztec diamond graph: Given an Aztec diamond graph of order 
n whose edges are marked with weights, first decompose the 
graph into 4-cycles (to be called "cells" following Mihai 
Ciucu's coinage [|l|l), and then replace each marked cell 




by the marked cell 

z/{wz + xy) X \^y/{wz + xy) 



x/ (wz + xy) 




w/{wz+xy) 



and strip away the outer flank of edges, so that an edge-marked 
Aztec diamond graph of order n — 1 remains. Then (as will be 
shown in section 5) the sum of the weights of the matchings 
of the graph of order n equals the sum of the weights of the 
matchings of the graph of order « — 1 times the product of the 
(different) factors of the form wz + xy ("cell-factors"). So, 
if you perform the reduction process n times, obtaining in the 
end the Aztec diamond of order (for which the sum of the 
weights of the matchings is 1), you will find that the weighted 
sum of the matchings of the original graph of order n is just 
the product of 



n2+(n-l)2 + . 



terms, read off from the cells. 



1' 
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Let us try this with the weighted Aztec diamond graph dis- 
cussed earlier, in connection with the 4x4 grid: 




First we apply the replacement rule described above: 




Next we strip away the outer layer (note that some of the val- 
ues computed at the previous step are simply thrown away, so 
that it is not in fact strictly necessary to have computed them): 




Now we do it again (only this time I will not bother to depict 
the weights that end up getting thrown away): 




We do it one more time, and we obtain the empty Aztec dia- 
mond. 

Now we multiply together all the cell-factors: 

2^ X (3/4)"* X (32/9) 36. 

A different example arises when the cells of an Aztec dia- 
mond of order n are alternately colored white and black, with 
all the edges in the white cells being assigned one weight (a) 
and all the edges in the black cells being assigned another 
weight (b). In this case, a single application of the algorithm 
gives an Aztec diamond of order n — 1 in which there are 
again edges of two different weights, only now the different- 
weighted edges are intermixed in the cells. However, if one 
repeats the algorithm again, one gets an Aztec diamond of or- 
der « — 2 in which the alternation of the two weights is as in 
the original graph. In this fashion, it is possible to express 
the sum of the weights of all the matchings of the graph in 
terms of cell-factors equal to 2a^, 2b^, and + b^. Indeed, 
by specializing these weights with a = 1 and b = we get 
a very simple proof of the "powers-of-5" formula proved by 
Yang [22 1, similar to the proof given by Ciucu [0]. 



3. COMPUTING EDGE-INCLUSION PROBABILITIES : 
THE ALGORITHM 

Let us assume that non-negative real weights have been as- 
signed to the edges of an Aztec diamond graph, and let us 
define the weight of a matching of the graph as the product of 
the weights of the constituent edges. Assume that the weight- 
ing of the edges is such that there exists at least one match- 
ing of non-zero weight. Then the weighting of the matchings 
induces a probability distribution on the set of matchings of 
the graph, in which the probability of a particular matching 
is proportional to the weight of that matching. This is most 
natural in the case where all edges have weight or weight 1 ; 
then one is looking at the uniform distribution on the set of all 
matchings of the subgraph consisting of the edges of weight 
1 . In any case, the probability distribution is well-defined, so 
it makes sense to ask, what is the probability that a random 
matching (chosen relative to the aforementioned probability 
distribution) will contain some particular edge e? Here we 
present a scheme that simultaneously answers this question 
for all edges e of the graph. 

What we do is thread our way back through the reduction 
process described in the previous section, starting from order 
and working our way back up to order n, computing the 
edge probability in successively larger and larger weighted 
Aztec diamonds graphs; along the way we make use of the 
cell-weights that were computed during the reduction pro- 
cess. Suppose we have computed edge probabilities for the 
weighted graph of order n — L To derive the edge probabili- 
ties for the weighted graph of order n, we embed the smaller 
graph in the larger (concentrically), divide the larger graph 
into cells, and swap the numbers belonging to edges on 
opposite sides of a cell (the example below will make this 
clearer). When we have done this, the numbers on the edges 
are (typically) not yet the true probabilities, but they can be 
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regarded as approximations to them, in the sense that we can 
add a slightly more complicated correction term that makes 
the "approximate" formulas exact. Let us zoom in on a partic- 
ular cell to see how it works: Consider the typical cell 



cell: 




and suppose that the edges with weights w, x, y, and z. (before 
the swapping has taken place) have been given "approximate" 
probabihties p, q, r, and s, respectively. Then the exact prob- 
abilities for these respective edges are 



^ {1 — p — q — r — s)wz ^ {I — p — q ~ r~ s)xy 



wz + xy 
(1 — p — q — r — s)xy 

q + - — ^^^7 — P + 

wz+xy 



wz+xy 
(1 — p~q — r — s)wz 
wz + xy 



The number \ — p — q — r— s will be called the deficit asso- 
ciated with that cell, and in the context of the four preceding 
inset expressions, it is called the (net) creation rate; the num- 
bers wz/{wz + xy) and xy/ {wz + xy) are called (net) creation 
biases. 

Applying this to the 4x4 grid graph, we find that for the 
weighted Aztec graphs of orders 1, 2, and 3 that we found 

(in reverse order) when we counted the matchings, the edge- 
probabilities are as follows: For the order 1 graph, we have 




We embed this in the Aztec graph of order 2: 




Note that I have put a dot in the middle of each of the four 
cells. We swap each number with the number opposite it in its 




These are the approximate, inexact edge-probabilities. To find 
the exact values, we compute that in each of the four cells the 
deficit I— p — q — r — sis, 1/2. The weights on the edges (all 
equal to 1 or 1 /2 — see second-to-Ust figure of section 2) give 
us creation biases 



and 



(l/2)/(3/4) = 2/3 



(l/4)/(3/4) = l/3. 



So we increment the 1 /2's (and the O's that are opposite them) 
by 

(l/2)(2/3) = l/3, 
and we increment the other O's by 

(l/2)(l/3) = l/6, 

obtaining 




Embed this in the graph of order 3, and swap the numbers in 
each cell: 
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The four corner-cells have deficit 1 /6, the other four outer 
cells have deficit 2/3 and the inner cell has deficit —1/3 (or 
surplus +1/3). So we do our adjustments, this time with all 
creation biases equal to 1/2 (except in the corners, where the 
creation biases are and 1 — see third-to-hst figure of section 
2): 




These are the true edge-inclusion probabilities associated with 
the original graph. That is, if we take a uniform random 
matching of the 4x4 grid, these are the probabilities with 
which we will see the respective edges occur in the matching. 

4. GENERATING A RANDOM MATCHING: THE 
ALGORITHM 

Before one can begin to generate random matchings of a 
weighted Aztec diamond graph of order «, one must apply the 
reduction algorithm used in the counting algorithm, obtain- 
ing weighted Aztec diamond graphs of orders n — 1, n — 2, 
etc. But once this has been done, generating random match- 
ings of the graph is quite easy. The algorithm is an iterative 
one: starting from a Aztec diamond graph of order 0, one 
successively generates random matchings of the graphs of or- 
ders 1, 2, 3, etc., using the weights that one computed during 
the reduction-process. At each stage, the probability of seeing 
any particular matching can be shown to be proportional to the 
weight of that matching (where the constant of proportionality 
naturally changes as one progresses to larger and larger Aztec 
diamond graphs). 

Here is how one iteration of the procedure works. Given a 
perfect matching of the Aztec diamond graph of order k — 1, 
embed the matching inside the Aztec diamond graph of order 
k, so that you have a partial matching of an Aztec diamond 
graph order k. (Note that the 4-cycles that were cells in the 
small graph become the holes between the cells in the large 
graph). In the new enlarged graph, find all (new) cells that 
contain two matched edges and delete both edges. (In [|||, this 
was called "destruction".) Then replace each edge in the re- 
sulting matching with the edge opposite it in its cell. (In [|^], 
this was called "shuffling", although in subsequent talks and 
articles I prefer to call this step of the algorithm "sliding" and 
to reserve the term "shuffling" for the compound operation 



consisting of "creation", "sliding", and "destruction".) It can 
be shown that the complement of the resulting matching (that 
is, the graph that remains when the matched vertices and all 
their incident edges are removed) can be covered by 4-cycles 
in a unique way (and it is easy to find this covering just by 
making a one-time scan through the graph). Each such 4-cycle 
has weights attached to each of its 4 edges, so we can choose 
a random perfect matching of the 4-cycle using the weights 
to determine the respective probabilities of the two different 
ways to match the cycle. (In this was called "creation".) 
Taking the union of these new edges with the edges that are 
already in place, we get a perfect matching of the Aztec dia- 
mond graph of order k. 

Let us see how this works with the 4x4 grid (for which we 
have already worked out the reduction process). To generate 
a random perfect matching, we start with the empty matching 
of the graph of order 0, embed it in the graph of order 1 , and 
find that the complementary graph of the empty matching is a 
single 4-cycle, which we match in one of two possible ways. 
In this case, each of the four edges has weight 4/3, so both 
matchings have equal likelihood; say we choose 




(Here we have used open circles to denote the vertices of the 
graph.) Now we embed this matching in the Aztec diamond 
graph of order 2: 

O O 




O O 

Neither of the two edges needs to be destroyed, because they 
belong to different cells. Replacing each by the opposite edge 
in its cell, we get 
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The complementary graph has a unique cover by 4-cycles: 



4-cycle, obtaining 





o 




o 



o 



o 



o 



o 



o 



o 



o 



hi each 4-cycle, three edges have weight 1/2 and one has 
weight 1, so we are twice as likely to see the one of the 
matchings as the other. Let us suppose that when we choose 
matchings with suitable bias, we get the more likely of the two 
matchings in both cells (as happens 4/9 of the time). Thus we 
have 




O 



o • o 



o 



o 




There is a unique cover of the complement by 4-cycles: 
O A O 





Now we embed the matching in the graph of order 3: 



o 



o 



o 



As it happens, each edge in the complement has weight 1 in 

the (original) weighted Aztec diamond of order 3, so we can 
use a fair coin to decide how to match each of the 4-cycles. 



O 



O 



O 




O 



• O 



O 



5. COMPUTING WEIGHT-SUMS: THE PROOF 

To prove the claim, we begin with a lemma (the "urban re- 
newal lemma"): If you have a weighted graph G that contains 
the local pattern 



o 



o 



o 



This time we must delete the two central edges, because they 
belong to the same cell. After we have done this, we replace 
each of the four remaining edges with the opposite edge of its 
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(here A,B,C,D are vertex labels and w,x,y,z are edge weights, 
with unmarked edges having weight 1) and you define G' as 
the graph you get when you replace this local pattern by 




with 

w' =z/{wz + xy) x' =y/{wz + xy) 

y' =x/{wz+xy) z! =w/{wz + xy), 

then the sum of the weights of the matchings of G equals 
wz + xy times the sum of the weights of the matchings of G' . 
The proof will be by verification: For each possible subset S of 
{A,B,C,Z)}, we will check that the total weight of the match- 
ings of G in which the vertices in 5 are matched with vertices 
inside the patch and the vertices in {A,fi, C,Z)} \5 are matched 
outside the patch equals wz+xy times the total weight of the 
matchings of G' in which the vertices in S are matched with 
vertices inside the patch and the vertices in {A , B, C,D} \ 5 are 
matched outside the patch. As it turns out, 10 of the 2^ = 16 
cases are trivial, and of the remaining 6, four are related by 
symmetry, so it comes down to three cases. 

• S = {A,B}: In this case our matching of G matches A,B 

inward and C,D outward, so it must contain the edge 
of weight z- This matching of G is associated with a 
matching of G' in which A and B are matched (inward) 
to each other via an edge of weight w' and C and D 
are matched outward just as before. The ratio of the 
weight of the (given) matching of G to the (associated) 
matching of G' is z/w' = wz + xy. 

• S= {A,B,C,D}: In this case our matching of G matches all 

four vertices A,B,C,D inward. This matching is asso- 
ciated with two matchings of G', one of which matches 
A with B and C with D and the other of which matches 
A with C and B with D (all edges outside of the patch 
remain as before). If we let Q denote the weight of the 
given matching of G, one of the two derived matchings 
of G' has weight w'z'Q and the other has weight a^/Q. 
The combined weight of the two matchings is 

{w'z!+xfy')Q = Q/ {wz+xy), 

so the ratio of the weight of the single (given) matching 
of G to the two (associated) matchings of G' is wz + xy. 



• 5 = 0: In this case our matching of G matches all four ver- 
tices A, B, C,Z) outward, and either contains edges w and 
z or contains edges x and y. These two cases must be 
lumped together. Let Q denote the product of the edges 
of all the weights of the other edges of the matching 
(not counting w,x,y,z). We find that two matchings of 
G, of weight wzQ and xyQ, are associated with a single 
matching of G' of weight Q. The ratio of the weights of 
the two matchings of G to the single matching of G' is 
wz + xy. 



This completes the proof of the lemma. 

Now, to prove the Aztec reduction theorem, suppose we 
have a weighted Aztec graph of order n. 




(Here n = 3 and the weights are not shown.) We begin our 
reduction by spUtting each vertex into three vertices: 




It is easy to see that this vertex- splitting does not change the 
sums of the weights of the matchings, as long as the new edges 
we add are given weight 1. 
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Now we apply urban renewal in locations to get a graph 
of the form 




with new weights w' , x', y' , z! replacing the old weights w, x, 
y, z- (If one were attempting a literal description of the algo- 
rithm, one would want the weight-variables to have subscripts 
i,i ranging from 1 to n.) Here is the coup de grace: The pen- 
dant edges that we see must belong to every perfect matching 
of the graph, so we can delete them from the graph, obtaining 
an Aztec graph of order n — 1 : 




The theorem now follows. 

One can prove the "fortress" theorem (see |^]) by similar 
means. Here, the initial graph is 




One's first impulse is to remove the pendant edges, but it is 
better still to apply urban renewal to every other "city" in 



this graph. Then one gets an Aztec diamond graph in which 
roughly half of the edges have weight 1 and the rest have 
weight 1 /2, with the two sorts of cities alternating in checker- 
board fashion. The same techniques that were used above can 
also be applied here to allow one to conclude that the number 
of matchings of such a graph is a certain power of 5 (or twice 
a certain power of 5). The "5" comes from the cell-factor 
(1)(1) + (l/2)(l/2) (when multiplied by 4), just as (in the 
unweighted Aztec diamond theorem) the "2" comes from the 
cell-factor (1)(1) + (1)(1). In short, the proof of Yang's theo- 
rem becomes a near-triviality (just as it is for Ciucu's method 
il). 

6. COMPUTING EDGE-INCLUSION PROBABILITIES : 
THE PROOF 

In this section I will prove that the iterative scheme for com- 
puting these probabilities that was described earlier does in 
fact work (at least when all edge-weights are non- vanishing). 
To start things off we will need a lemma: This lemma is a gen- 
eralization of a proposition originally conjectured by Alexan- 
dru lonescu in the context of Aztec diamond graphs, and was 
proved by Propp in 1993 (private e-mail communication). 

Let G be a bipartite planar graph and let A,B,C,D be four 
vertices of G that form a 4-cycle ABDC bounding a face of G: 




Assume the edges of this graph are weighted in some fashion. 
Define G^g as the weighted graph obtained by deleting ver- 
tices A and B and all incident edges; all other edges have the 
same weight as in G. Define G^^, Gg^, G^^, and G^^^^ sim- 
ilarly. For any edge-weighted graph H, let M{H) denote the 
sum of the weights of the matchings of H. Then the identity 
we will prove is 

M{G)M{G^^cd) ^M{G^b)M{Gco) +M{G^c)M{Gso). 

Note that if we superimpose a matching of G and a match- 
ing of Gj^gfjij, we get a multiset of edges of G in which ver- 
tices A,B,C,D get degree 1 and every other vertex gets degree 
2. The same thing happens if we superimpose a matching 
of G^g and a matching of G^^, or a matching of G^^ and a 
matching of G^^,. Call an edge-multiset of this sort a "near 2- 
factor" of G. I will prove that equality holds above by showing 
that each near 2-factor F makes the same contribution to the 
left side as it does to the right side. 
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Pick a specific near 2-factor F of G; it consists of k 
closed loops plus some paths that connect up the four ver- 
tices A,B,C,D in pairs. It is impossible that the vertices are 
paired as {A,D} and {B,C}, since the graph is planar. Hence 
the vertices pair as either {A,B} and {C,D} or {A,C} and 
{B,D}. Without loss of generality, let us focus on the former 
case. Since G is bipartite, the path from A to B will contain an 
odd number of edges, as will the path from C to D. It is easy 
to verify the following three claims: 

1 . There are 2*^ ways to split up F as a matching of G to- 
gether with a matching of G^gc^. 

2. There are 2*^ ways to split up F as a matching of G^g 
together with a matching of G(^jj. 

3. There are no ways to split up F as a matching of G^^ 
together with a matching of Ggj^. 

Moreover, one can check that each near 2-factor that pairs the 
vertices as {A,B} and {C,D} contributes equal weight to both 
sides of the equation. Clearly the same is true for the near 
2-factors that pair the vertices the other way (A with C and B 
with D). This completes the proof. 

(Note: Another way to prove this lemma is to use the face 
that the number of perfect matchings of a bipartite plane graph 
with 2« vertices can be written as the determinant of an n-by-« 
Kasteleyn-Percus matrix [ [l2t [[ist. Then the lemma is seen to 
be a restatement of Dodgson condensation see also jl^). 

It will be convenient for us to restate this lemma in a some- 
what different form. Assume that the graph has at least one 
matching of positive weight. Let w, x, y, and z denote the re- 
spective weights of the edges AB, AC, BD, and CD, and let p, 
q, r, and s denote the probabilities of these respective edges 
being included in a random matching of G (where as usual 
the probability of an individual matching is proportional to its 
weight): 




probability that a random matching of G contains both edge 
AB and edge CD is {psxy + qrwz) /xy, while the correspond- 
ing probabiUty for edges AC and BD is {psxy + qrwz) / wz. (To 
see that this follows from the formula for P*, note that their 
sum is indeed {wz + xy){ps / wz + qr / xy) — P* , and that their 
ratio is wz/xy, which is the ratio of the weights of any two 
matchings that differ only in that one contains AB and CD 
while the other contains AC and BD.) 

Now let us go back and give urban renewal a fresh look. 
Here are the graphs G and G', with edges marked with both 
their weight and their probability (in the format "weight [prob- 
ability]"). Here is G: 




(The unmarked edges have weight 1, and their probabilities 

are I — p — q, 1 — ^ — et cetera.) 
Here is G': 




y [r]\ / z [s] 

DX/ 



The probability P* that a random matching of G has 
an alternating cycle at face ABDC (i.e., either includes 
edges AB and CD or includes edges AC and BD) is 
equal to M{G^gi^ij){wz + xy)/M{G). Since we have p = 
wM{G^„)/M{G), q=xM{G^c)/M{G),r^yM{Gso)/M{G), 
and s — zM{G(ju)/M{G), the preceding lemma tells us that, 
in the event that w,x,y,z are all non-zero, the probability P* is 
equal to {psxy + qrwz){-^ + It)' follows from this that the 



Recall that we are in the situation where we know p', q', r', 
and s', and are trying to compute p, q, r, and s. Recall also 
that 

w' —z/{wz + xy) x! —y/{wz + xy) 

y' —xl{wz + xy) z! = w/ {wz+xy). 

Here is our plan: First, we will work out the probabilities 
of all the local patterns in G'. Then we will use the urban re- 
newal correspondence to deduce the probabilities of the local 
patterns in G. Finally, we will deduce the probabilities of the 
individual edges in G. 

1 will write p' ,q' ,r' ,s' ,w' ,x' ,y' ,z! as P,Q,R,S,W,X,Y,Z, to save 
eye-strain. 1 will also write A' = XYPS + WZQR. 
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Let S denote the set of vertices in {A,B,C,D} that are 
matched to a vertex inside the patch. Then here are the re- 
spective probabilities of the local patterns in G': 

S^{A,B,C,D} : {WZ + XY)A'/WXYZ 
S = {A,B}:P-A'/XY 
S^{A,C}:Q-A'/WZ 
S = {B,D}:R-A'/WZ 
S = {C,D}:S~A'/XY 

S = {}: l-P-Q-R-S+{WZ + XY)A'/WXYZ 

(The second formula is obtained by recalling that the proba- 
bility of edge AB being present in a random matching, which 
can happen in two ways according to whether or not edge CD 
is present, must be P; the last formula is obtained by recalling 
that the probabilities of the six local configurations must sum 
to 1.) 

It follows from the last of these, and the urban renewal cor- 
respondence, that the probability that a random matching of G 
contains two edges from the 4-cycle ABCD is 

l-P-Q-R-S+{WZ + XY){PS/WZ + QR/XY). 

The probability that a random matching of G contains edges 
AB and CD must be equal to this quantity times 



wz/(wz + xy) ^ WZ/{WZ + XY), 



which gives 



{WZ/{WZ + XY)){\-P-Q-R-S)+PS + QR{WZlXY). 

The probability that a random matching of G contains edge 
AB but not edge CD (by another application of the urban re- 
newal correspondence) equals the probability that a random 
matching of G' contains edge CD but not edge AB, which is 

S-PS-QR{WZ/XY). 

Adding, we find that the probability that a random matching 
of G contains edge AB is 



wz 



wz+xy 
Thus we conclude that 

wz 



wz + xy 



(l-P-Q-R-S). 



(l-P-Q-R-S), 



where wz/ (wz+xy) is the creation bias and I— P Q R S 
is the deficit. The formulas for q, r, and s follow by symmetry. 



7. GENERATING A RANDOM MATCHING: THE PROOF 

Now I will show that (generalized) domino-shuffling 
works. 

First, let us note that the algorithm for generating random 
matchings that was presented earlier really is a generalization 



of the domino-shuffling procedure described in [||]. The op- 
eration of removing from a matching those matched edges 
that belong to the same cell as another matched edge corre- 
sponds to the process of removing odd blocks to obtain an 
odd-deficient tiling; the operation of swapping the remaining 
edges to the opposite side of the cell corresponds to the pro- 
cess of shuffling dominos (resulting in the creation of an odd- 
deficient tiling of a larger Aztec diamond); and the operation 
of introducing new edges corresponds to the creation of new 
2x2 blocks. 

Next, let us recall that a perfect matching of the Aztec di- 
amond graph of order n determines an n x « alternating-sign 
matrix, as was first described in [|^]. Specifically, for each 
of the cells of the graph, write down the number of edges 
from that cell that participate in the matching and subtract 1 . 
For instance: 




Label the weights of the edges of the cell in row i and column 
/ with W;, ■, X ;, V ;, and z- ., thus: 




Define also a "cell-weight" 

D,j = ^,fy+Xijyij 

(these weights are the same as the cell-factors considered ear- 
lier in this article). Call two edges in a matching "parallel" if 
they belong to the same cell. 

The bias in the creation rates for generalized shuffling was 
chosen so that the ratio between the probabilities of choos- 
ing one pair of parallel edges versus another (wz versus xy) 
is equal to the ratio of the weights of the resulting matchings. 
Consequently, if two perfect matchings of the Aztec graph of 
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order n are associated with the same « x n alternating-sign ma- 
trix, then their relative probabilities are correct. So all we need 
to do is verify that the aggregate probability of these match- 
ings (for fixed alternating-sign matrix A) is what it should be. 

To do this, it is convenient to introduce certain partial 
matchings of Aztec diamond graphs (analogous to the odd- 
deficient and even-deficient domino-tilings considered in [|[]). 
For the Aztec diamond of order «, we consider partial match- 
ings obtained from perfect matchings by removal of all par- 
allel edges; that is, pairs of edges that belong to the same 
cell. For the Aztec diamond of order n — 1, we consider par- 
tial matchings obtained from perfect matchings by removal of 
all pairs of edges that belong to the same "co-cell" (defined 
as the space between four adjoining cells). The "destruction" 
phase of generalized shuffling converts a perfect matching of 
the graph of order n — 1 into a partial matching M of the graph 
of order « — 1, the "sliding" phase converts that partial match- 
ing into a partial matching M' of the graph of order n, and 
the "creation" phase converts that partial matching into a per- 
fect matching. Note that M determines M' uniquely, and vice 
versa. 

Define the aggregate cell-weights D^, Dg, and D to be 
the products of those cell-weights j associated with loca- 
tions in the Aztec diamond at which +1,0, and — 1 appear, re- 
spectively (so that D^DqD is just the product of all the cell- 
weights). The following three claims can be readily checked: 
the sum of the weights of the perfect matchings of the Aztec 
diamond of order n that extend the partial matching M' is 
equal to times the weight of M' itself (defined as the prod- 
uct of the weights of its constituent edges); the weight of M' 
is equal to Dq times the weight of M (defined in terms of the 
weighting of the graph of order n — 1); and the weight of M is 
equal to D times the sum of the weights of the perfect match- 
ings of the Aztec diamond of order « — 1 that extend the partial 
matching M. (This relationship between cell-weights and the 
entries of alternating-sign matrices was originally made clear 
in the work of Ciucu [|l|].) 

The upshot is, the sum of the weights of the perfect match- 
ings of the Aztec diamond graph of order n that are associated 
with the alternating- sign matrix A equals 

^1.1^1 2 ■ ■ '^"J^ 

times the sum of the weights of all the perfect matchings of the 
Aztec diamond graph of order n — 1 that are capable of giving 
rise to it under shuffling. This factor is independent of A. It 
follows that if one takes a random matching of the smaller 
graph (in accordance with the edge weights given by urban 
renewal) and applies destruction, shuffling, and creation, then 
one will get a random matching of the larger graph. 

8. DIABOLO-TILINGS OF FORTRESSES 

Consider an n-by-n array of unit squares, each of which 
has been cut by both of its diagonals, so that there are 
isosceles right triangles in all. A region called a fortress of 
order n is obtained by removing some of the triangles that are 
at the boundary of the array. More specifically, one colors 



the triangles alternately black and white, and removes all the 
black (resp. white) triangles that have an edge on the top or 
bottom (resp. left or right) boundaries of the array. Here, for 
instance, is a fortress of order 4: 




(When n is even, the two ways of coloring the triangles lead 
to fortresses that are mirror images of one another; however, 
when n is odd, the two sorts of fortresses one obtains by re- 
moving the colored triangles in the specified fashion are gen- 
uinely different.) The graph dual to a fortress (with vertices 
corresponding to the isosceles right triangles) is composed of 
4-cycles and 8-cycles. (If one removes the pendant edges from 
the final figure in section 5 one obtains the graph dual to a 
fortress of order 3.) 

We call the small isosceles right triangles monobolos. 
Given our tiling of the fortress by monobolos, we can (in 
many different ways) form a new tiling with half as many 
tiles by amalgamating the monobolos in pairs, where the two 
monobolos that are paired together are required to be adja- 
cent. These new tiles are diabolos, and come in two different 
shapes: squares and isosceles triangles. (Note: in the liter- 
ature on tilings, such as [^], a third kind of diabolo is rec- 
ognized, namely a parallelogram with angles of 45 and 135 
degrees with side-lengths in the ratio 1 : \/2; however, these 
diabolos do not arise in the context being discussed here.) 

A tiling obtained from the original monobolo-tiling of the 
fortress by amalgamating adjacent pairs of monobolos will be 
called, for (comparative) brevity, a diabolo-tiling of a fortress. 
It should, however, be kept in mind that we are tacitly limiting 
ourselves to those tilings that can be obtained by the afore- 
mentioned process of amalgamation; or, if one prefers to use 
the dual picture, tilings that correspond to perfect matchings 
of the "squares and octagons" graph that is dual to the original 
tiling by monobolos. For instance, the first of the following 
two tilings of a fortress of order 3 by means of tiles that are 
diabolos is an example of what we mean by a "diabolo-tiling 
of a fortress", but the second is not, because the four center 
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diabolos are not amalgams of the required kind: 




To count the diabolo-tilings of a fortress (understood in the 
above sense), one replaces the underlying monobolo-tiling by 
its dual graph and adds pendant edges along the boundary, 
obtaining a graph like the one shown at the end of section 
5. One can then apply the urban renewal lemma to turn this 
into an Aztec diamond graph with edges of weight 1 arranged 
in cells that alternate, checkerboard-style, with cells whose 
edges have weight j. In this fashion one can derive Yang's 
formula. 

Now, however, we set our sights higher. In section 6, we 
gave a scheme for calculating the probability that a given 
edge of a weighted Aztec diamond graph is included in a ran- 
dom matching, where the probability of a matching is propor- 
tional to its weight. This in fact lets us write down an explicit 
power series whose coefficients are precisely those inclusion- 
probabilities. We will show how this goes in a slightly more 
general setting, where the edges of weight j are replaced by 
edges of weight f . 

Here we will find it convenient to use a rotated and shifted 
version of the Aztec diamond graph of order n, whose edges 
are all either horizontal or vertical (rather than diagonal). The 
vertices of this graph are the integer points where |/ + 
2 1 + U + 7 1 ^ the edges of the graph join vertices at distance 
1. The horizontal edge joining (/ — and is called 



"north-going" if / + j + n is odd and "south-going" otherwise, 
and the vertical edge joining 1) and is called "east- 
going" if / + j + n is odd and "west-going" otherwise (this 
nomenclature is taken from [|||). The kind of weighting we 
are considering can be described as follows: 

1 . Each horizontal edge has weight 1 or weight t accord- 
ing to the parity of its x-coordinate. (That is, if two 
such edges are related by a unit vertical displacement, 
they have same weight, but if they are related by a unit 
horizontal displacement, their weights differ.) 

2. Each vertical edge has weight 1 or weight t according 
to the parity of its y-coordinate. 

3. If n is 1 or 2 (mod 4), the four extreme-most edges (the 
northernmost and southernmost horizontal edges and 
the easternmost and westernmost vertical edges) have 
weight f ; if n is or 3 (mod 4), those four edges have 
weight 1. 

(These four situations are considered together because each 
is linked to the next by the shuffling operation. Technically, 
this scheme only handles fortresses of odd order, and indeed, 
only handles half of those, but the other case are quite similar) 
We define a generating function P{x,y,z) in which the coeffi- 
cient of jc'yz" (with k > 0, i + j + k odd, and |/| + \j\ < n) is 
the probabihty that a random matching of the Aztec diamond 
graph of order n, chosen in accordance with the probability 
distribution determined by the weight t, contains a horizon- 
tal bond joining vertices (/— and where / + j + n 
is odd. (We are limiting ourselves to northgoing horizontal 
bonds but this limitation is unimportant, since by symmetry 
the other three sorts of bonds have the same behavior) Note 
that the coefficients in this generating function (viewed as a 
power series in the variables x, y, and z) are polynomials in f , 
and that for each n, only finitely many pairs /,y contribute a 
non-zero term to P{x^y,z). 

The presence of t in the pattern of edge-weighting makes it 
natural to divide the north-going edges into four classes, ac- 
cording to how they "sit" relative to the weighting (cells are 
weighted in two patterns when n is odd and in two other pat- 
terns when n is even). Correspondingly, we write P{x,y,z) 
as a sum of four generating functions, each corresponding to 
one of the four classes of north-going edges. We call these 
generating functions E, F, G, and H. We also introduce gen- 
erating functions whose coefficients are the net creation rates 
of the four kinds of differently-weighted cells, which we de- 
note by e, /, g, and h. E and e are associated with cells in 
which all edges have weight f; F and / are associated with 
cells in which all edges have weight 1 ; G and g are associated 
with cells in which horizontal edges have weight t and vertical 
edges have weight 1 ; and H and h are associated with cells in 
which horizontal edges have weight 1 and vertical edges have 
weight t. Generalized shuffling implies the following alge- 
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braic relationships among the four generating functions: 



i + T-z7^ii^+^'')8 + iy+y'')h)-z^f 



i+t 



;z{{x + x ')h+{y + y ^)8)-z^" 



l+t- 



-z{{x + x-')f+{y + y-')e)~z^h 



1 



h^-z{{x + x-y+{y + y-')f)~z\ 

These can be solved simultaneously, giving us four three- 
variable generating functions e{x,y,z), f{x,y,z), g{x,y,z), 
h{x,y,z). We can then solve 

E^Hyz+^ez F^Gyz+^fz 

f2 1 

G = Eyz + ^ , ,9 H = Fyz+-r—^hz 



l+f2' 



and obtain (messy!) generating functions E{x,y,z), F{x,y,z), 
G{x,y,z), H{x,y,z)- The sum E +F + G + H, multiplied by z, 
is the desired generating function P. It is too messy to write 
down here, but it does at least in principle give one a way to 
understand what is going on with random diabolo-tilings of 
large fortresses. 

We can also use generalized shuffling to generate a random 
tiling of a large fortress. Recall that when this was done for 
Aztec diamonds in the early 1990s, a new phenomenon was 
discovered: the "arctic circle" effect jlTI]. Specifically, with 
probability going to 1 as the size of the Aztec diamond graph 
increases, a random tiling divides itself up naturally into five 
regions: four outer, "frozen" regions in which all the domi- 
nos line up with their neighbors in a repeating pattern, and a 
central, "temperate" region in which the dominos are jumbled 
together in a random-looking way. One feature of the tem- 
perate zone, rigorously proved in [|]|, is that the asymptotic 
local statistics, expressed as a function of normalized posi- 
tion, are nowhere constant; that is, the probability of seeing a 
particular local pattern of dominos in a random tiling changes 
as one shifts one view over macroscopic distances (i.e., dis- 
tances comparable to the size of the Aztec diamond graph). 
Strikingly, the boundary of the temperate zone, rescaled, tends 
in probability to a perfect circle. 

One might expect random diabolo-tilings of fortresses to 
exhibit much the same sort of phenomena, albeit with the cir- 
cular temperate zone probably replaced by a temperate zone 
of some other shape. However, when the generalized shuf- 
fling algorithm was used to generate random tilings of large 
fortresses, a startling new phenomenon appeared: within the 
very inmost part of the fortress, local statistics do not appear 
to undergo variation. That is to say, within this region, dubbed 
the "tropical astroid", the local statistics of a random tiling ap- 
pear to be shielded from the boundary, so that the (normalized) 
position of a tile relative to the boundary does not matter (as 
longas it stays fairly close to the center of the region). Fig- 
ure ll shows a random diabolo-tiling of the fortress of order 
200; the square diabolos have been shaded, to highlight the 



shape of the curve that separates the tropical zone from the 
temperate zone. 

Suppose we put coordinates on the fortress so that the cor- 
ners are at (0,2), (0,-2), (2,0), and (—2,0). Then it appears 
that the boundary of the frozen region is given by one real 
component of the curve 

400jc^ + 400y^ + 3400x2y'^ + 3400y^x'' + 8025jcy 
+ IOOOjc*' + lOOOy*^ - lllSOxY - 17250xy - 143 Ix^ 
- 143 1/ + 258 12xy - 3402x2 - 3402/ + 729 = 

(the "octic circle") and that the boundary of the tropical region 
is given by the other real component. Henry Cohn and Robin 
Pemantle, as of this writing, are working on a proof of these 
two assertions, by giving an asymptotic analysis of the coef- 
ficients of the generating functions obtained via generalized 
shuffling. 



9. LAST THOUGHTS 

Here are some last thoughts concerning urban renewal, 
domino-shuffling et cetera. 

1. This article has glossed over an important point, namely, 

anomalies that arise when one attempts to apply the al- 
gorithm to Aztec diamond graphs in which some edges 
have been assigned weight 0. Even if the weighted 
Aztec diamond in question has matchings of positive 
weight, it is possible that somewhere in the reduction 
process one will encounter cells with cell- weight 0; this 
leads to blow-up problems when one attempts to divide 
by the cell-weight. In such a case, one should replace 
edge- weights equal to in the original graph with edge- 
weights equal to e; the result will be a rational function 
of e whose limiting value as e goes to zero is desired. 
For this purpose, every rational function of e can be 
written as a polynomial or power series in e and re- 
placed by its leading term, so calculations are not as 
bad as one might think. This will work as long as the 
original weighted Aztec diamond graph has at least one 
matching of positive weight. 

2. These algorithms arose partly in response to work of Ku- 

perberg (see |16|), which took a more algebraic per- 
spective on enumeration, using the approach pioneered 
by Kasteleyn [|l2|. Kasteleyn's method requires mak- 
ing some arbitrary choices that end up not affecting 
the final answers to meaningful enumerative questions. 
This new work arose out of an attempt to find the in- 
variant combinatorial core of the Kasteleyn method. In 
particular, graph-rewriting is a combinatorial substitute 
for algebraic operations like row-reduction applied to 
a Kasteleyn matrix. However, this analogy was never 
worked out in any kind of rigorous detail. It would be 
helpful to see this explained. 



3. Continuing the above remark: One strength of the alge- 
braic approach, exploited by Kenyon [O] and others, is 
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that edge-inclusion probabilities, and, more generally, 
"pattern-occurrence probabilities", can be expressed in 
terms of determinants of minors of the inverse of the 
Kasteleyn matrix. Here again, the answer must be in- 
dependent of the arbitrary choices that were made in 
forming the matrix. So it is natural to hope that there 
will be a similar canonical scheme for calculating such 
pattern-occurrence probabilities as well. This may be 
similar to the problem of finding an extension of Dodg- 



son's condensation scheme that permits one to ef- 
ficiently compute the inverse of a matrix and not just 
compute its determinant. 

Generalized shuffling bears a strong resemblance to an al- 
gorithm described in recent work of Viennot |20|. I 
have not studied the matter deeply enough to be able 
to identify the relationship, but 1 strongly believe that 
the resemblance is more than coincidental. 
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5. It would also be desirable to extend urban renewal to 

matchings of non-bipartite planar graphs. This is equiv- 
alent to asking for routinized matrix-reduction schemes 
for Kasteleyn's Pfaffians. 

6. As described in [||, lonescu's recurrence gives rise to edge- 

inclusion probabilities for uniformly-weighted perfect 
matchings of Aztec diamond graphs, and it is observed 
that if two edges e, e' are close to one another (relative 
to the overall dimensions of the large Aztec diamond 
graph they both inhabit), and e and e' occupy identi- 
cal positions within their respective cells (i.e., both are 
either the northwest, northeast, southwest, or southeast 
edges in their cells), then the edge-inclusion probabili- 
ties for e and e' are very close numerically. Indeed, this 
phenomenon is a rigorously-proved consequence of the 
detailed analysis given in [^j]. However, it would be 
good to have a conceptual explanation for this continu- 
ity property. In particular, it seems plausible that the 
recurrence relation for edge-inclusion probabilities has 
the property of smoothing out differences. A clear for- 
mulation of such a smoothing-out property, and a rigor- 
ous proof that it holds, would be very desirable, since it 



might lead to a proof that this kind of continuity prop- 
erty holds for lozenge tilings of hexagons. (Numeri- 
cal evidence supports this continuity conjecture, but the 
method used in 1 1] does not permit one to draw conclu- 
sions of this nature.) On the other hand, it should be 
noted that these issues are subtle; for instance, net cre- 
ation rates (viewed as a function of position) tend not to 
vary as smoothly as edge-inclusion probabilities. 

7. As is explained in enumeration of perfect matchings of 
Aztec diamond graphs corresponds to "2-enumeration" 
of alternating-sign matrices, where the weight assigned 
to a particular ASM is 2 to the power of the number of 
— I's it contains. One can more generally consider x- 
enumeration, where 2 is replaced by a general quantity 
x; ordinary enumeration corresponds to 1 -enumeration, 
and 3-enumeration of ASMs leads to interesting exact 
formulas (see e.g. |]l7|]). Could shuffling be extended to 
give a scheme for sampling from the set of ASMs with 
uniform distribution or more generally an x-weighted 
distribution? 



[1] Mihai Ciucu, Perfect matchings of cellular graphs, / Alg. Corn- 
bin. 5 {1996), 87-103. 

[2] Mihai Ciucu, A complementation theorem for perfect match- 
ings of graphs having a cellular completion, /. Combin. Theory 
Ser. A81 (1998), 34-68. 

[3] Henry Cohn, Noam Elkies and James Propp, Local statistics for 
random domino ti lings of the Aztec diam ond, Duke Math. J. 85 
(1996), 117-166; Imath. CO/0 00 82 4 3 



[4] Henry Cohn, Michael Larsen and James Propp, The shape of 
a typical boxed plane partition, A^ew York J. Math. 4 (1998), 



137-165; 



math. CO/9801059 



[5] Charles J. Colboum, J. Scott Provan and Dirk Vertigan, A new 
approach to solving three combinatorial enumeration problems 
on planar graphs. Disc. Appl. Math. 60 (1995), 119-129. 

[6] S.J. Cyvin and I. Gutman, Kekule Structures in Benzenoid Hy- 
drocarbons (Lecture Notes in Chemistry 46), Springer- Verlag, 
1988. 

[7] C.L. Dodgson, Condensation of determinants, Proc. Royal Soc. 

London 15 (1866), 150-155. 
[8] Noam Elkies, Greg Kuperberg, Michael Larsen and James 

Propp, Alternating sign matrices and domino tilings J. Alg. 

Combin.! (1992), 111-132, 219-234; currently available from 



http : / / www . math . wise . edu/ ~propp/ 



[9] Martin Gardner, Polyhexes and polyabolos, chapter 11 in: 
Mathematical Magic Show, Vintage Books, 1978. 

[10] D. Grensing, L Carlsen and H.C Zapp, Some exact results for 
the dimer problem on plane lattices with non-standard bound- 
aries, Phil. Mag. A 41 (1980), lll-lSl. 

[11] William Jockusch, James Propp and Peter Shor, Random 
domino tilings and th e arctic circle theorem, preprint 1995; 



Mth.co/gsoioei 



[12] P.W. Kasteleyn, Graph theory and crystal physics, in: Graph 
Theory and Theoretical Physics, F. Harary, ed.. Academic 
Press, 1967. 



[13] Richard Kenyon, Local statistics of lattice 
dimers, Ann. Ins. H. Poincare Probab. Statist. 
33 (1997), 591-618; currently availa ble from 



http : / /topo . math . u-psud . f r /~kenyon/ 



[14] Richard Kenyon, The planar dimer model with boundary: a sur- 
vey, Directions in mathematical quasicrystals, 307-328, CRM 
Mon ogr Ser 13, Amer Math. Soc. 2000; currently avail able 



from :ittp : / / topo . math . u-psud. fr/ ~kenyon/. 
[15] Eric Kuo, personal communication; pre- 
liminary write-up currently avail able from 



http : / / www . OS . berkeley .edu/ --^ekuo/ 



[16] Greg Ku perberg, An exploration of the permanent-determinant 



method; nath . CO/ 98 10 91 



[17] Greg Kuperberg, Symmetry class of alternating-sign ma- 
trices under one roof, to appear in Ann. of Math.; 



Tiath. CO/0008184 



[18] Jerome Percus, One more technique for the dimer problem, / 
Math. Phys. 10 (1969), 1881-1888. 

[19] David P. Robbins and Howard Rumsey, Jr., Determinants and 
alternating sign matrices. Adv. in Math. 62 (1986), 169-184. 

[20] X.G. Viennot, A combinatorial interpretation of the quotient- 
difference algorithm, in: Formal Power Series and Algebraic 
Combinatorics, 12th International Conference, D. Krob, A.A. 
Mikhalev and A.V. Mikhalev, eds., Springer, 2000. 

[21] David Wilson, Determinant algorithms for random 
planar structures, in: Proceedings of the Eighth An- 
nual ACM-SIAM Symposium on Discrete Algorithms 
(SODA 1997), 258-267: currentlv available from 



http : / /www .research . microsoft . com/ ~dbwilson/ 
[22] Bo- Yin Yang, Three enumeration problems concerning Aztec 
diamonds, Ph.D. thesis. Department of Mathematics, Mas- 
sachusetts Institute of Technology, Cambridge, Massachusetts, 
1991. 

[23] C. Zeng, P.L. Leath and T. Hwa, Thermodynamics of Meso- 



18 



scopic Vortex Systems in 1 + 1 Dimensions, Phys. Rev. Letts. 
83 (1999), 4860^683. 



